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During  the  past  year  our  work  supported  by  the  ONR 
Contract  N00014-84-C-0353  achieved  many  significant  results 
concerning  ionospheric  generation  of  ELF/VLF  as  well  as  (JLF 
waves.  A  major  breakthrough  was  our  participation  in  an 

experiment  using  the  Arecibo  facility  in  which  a  preliminary 
proof  of  principle  of  two  of  the  relevant  concepts  was 
produced.  We  highlight  below  some  of  our  key  results  and 
refer  the  interested  reader  to  the  appendices  for  details."^ 
(i)  A  novel  process  for  generating  ELF/ULF  waves  by 

using  ionospheric  heating  was  discovered  and  its 
applicability  to  auroral  geometry  was  exmined.  The  process 
relies  on  modulated  E-region  heating  and  does  not  depend^^n 
the  preence  of  cunbient  ionospheric  currents. y  The  ELF/ULF 
signals  are  caused  by  oscillatory  ionospheric  currents 
spontaneously  generated  in  the  HF  heated  volume  by  the 
coupling  of  the  hot  spot  temperature  gradient  (VT)  to  the 
amibient  ionospheric  density  gradient  (vnxVT) .  The  mechanism 
can  account  for  the  remarkably  large  ULF  signals  observed  by 

the  Tromso  HF  modifications  facility.  The  details  can  be 

found  in  a  paper  published  in  Geophysical  Research  Letters, 
12,  279,  (1985)  and  Included  as  Appendix  I. 


(ii)  A  comprehensive  analytic  and  computational  study  of 
the  ELF  generation  by  stimulated  beat  excitation  in  a 
dissipative  medium  was  completed.  The  power  thresholds  that 
lead  to  complete  pump  decay  for  a  transmitter  with  30° 
incidence  at  various  ionospheric  heights  are  given  in  Table 
I  of  Appendix  II  along  with  the  required  pump  frequency  and 
the  resulting  low  frequency.  The  details  can  be  found  in  a 
paper  accepted  for  publication  in  the  Journal  of  Geophysical 
Research  and  included  here  as  Appendix  II. 

(iii)  We  participated  in  the  experimental  program  of  Dr. 
Ganguly  (Rice  University)  using  the  Arecibo  facility.  The 
experiment  produced  signals  at  3  and  5  Hz  of  the  order  of  1 
my  when  operated  under  conditions  appropriate  for  stimulated 
beat  excitation.  ULF  waves  with  comparable  amplitude  were 
also  produced  by  using  the  technique  mentioned  in  (i)  above. 
The  interaction  occurred  in  the  F-region  and  therefore  it 
was  different  than  the  current  modulation  schemes.  The  work 
will  be  presented  as  an  invited  paper  at  the  APS  Plasma 
Physics  Meeting  in  San  Diego,  November  1985.  The  details 
can  be  found  in  a  paper  submitted  to  Phys.  Rev.  Lett,  and 
included  as  Appendix  III. 

(iv)  The  structure  and  properties  of  the  ELF  plasma 
eigenmodes  in  the  polar  E-region  of  the  ionosphere  was 
examined  in  detail.  It  was  noted  that  for  altitudes  below 
150km  the  ion  neutral  collision  frequency  (vi)  is  larger 


than  the  ion  cyclotron  frequency  (Di)  thereby  preventing  the 
existence  of  proper  Alfven  or  magnetosonic  eigenmodes.  Our 
work,  however,  indicates  that  the  strong  ion  collisionality 
(Vi  >  fli)  freezes  the  ion  motion.  This  result  shows  that 
the  whistler  mode  becomes  a  helicon  right  hand  polarized 
mode  for  frequencies  WLp  <  vt-  This  is  similar  in  nature  to 
the  well  known  in  solid  state  physics  helicon  mode,  in  which 
case  the  lattice  rather  than  the  collisionality  freezes  the 
ion  motion.  Since  both  the  downconversion  efficiency  and 
the  coupling  of  the  ELF  signals  to  the  earth  ionosphere 
waveguide  favor  low  altitudes,  it  was  concluded  that  the 
nonlinear  coupling  should  be  to  the  helicon  rather  than  the 
Alfven  or  magnetosonic  mode.  A  draft  of  the  paper,  is 
included  as  Appendix  IV. 

(v)  Our  most  recent  work,  left  unfinished  due  to  the 
termination  of  the  contract,  focused  on  the  possibility  of 
generating  oscillatory  magnetic  moments  that  radiate  at 
ELF/VLF  frequencies  in  the  ionosphere  by  utilizing  resonance 
absorption  generated  magnetic  fields.  The  field  generation 
by  this  process  has  been  studied  both  theoretically 
(Bezzerides  et  al.,  Phys.  Rev.  Lett,  38,  663,  1977)  and 
experimentally  (DiVergilio  et  al.,  Phys.  Rev.  Lett.  3^,  541, 
1977;  Obenschain  and  Luhmann,  Jr.,  Phys.  Rev.  Lett.  311, 
1979)  for  the  case  of  laser  or  microwaves  interacting  with 
the  plasma.  The  saturated  value  of  the  self  generated  field 


by  radiowaves  with  frequency  f  =  Uq/2v  and  quiver  velocity 
V  =  eE/inWo  is  (Bezzerides  et  al.,  Phys.  Rev.  Lett.  38,  495, 
1977) 


corresponding  to  a  power  PvLF  the  modulating  frequency  f 
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The  magnetic  field  B  on  the  ground  at  a  radius  distance  r 
from  the  modulated  region  will  be 
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which  is  many  tens  of  db  above  noise.  Our  preliminary 
results  indicate  that  this  is  an  extremely  important  and 
efficient  process  and  should  be  studied  in  appropriate 
depth,  both  for  proof  of  principle  experiments  in  Alaska  as 
well  as  for  practical  Navy  communications. 
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Abstract.  It  Is  shown  that  anplltude 
Modulated  HP  heatinc  in  the  ionospheric  E- region, 
can  produce  substantial  power  in  the  ELP/ULP 
range  even  in  the  absence  of  asbient  ionospheric 
currents.  The  signals  are  associated  with 
spontaneous  generation  of  eagnetic  fields  caused 
by  the  coupling  of  the  hot  spot  teaperature 
gradient  to  the  aabient  ionospheric  density 
gradient  (Vn  x  VT).  The  eechanise  produces 
values  for  PCS  eicropulsations  consistent  with 
those  observed  by  Stubbe  and  KopAa  ( 1981 ) . 

Introduction 

Excitation  of  vlp/elf/UIP  waves  by  the 
interaction  of  powerful  HP  waves  in  the 
ionosphere  has  been  a  very  active  theoretical  and 
ej^rlaental  topic.  Up  to  now  the  Most 
successful  excitation  techniques  have  been 
Modulation  of  the  polar  and  equatorial 
electrojets.  In  this  technique  an  HP  wave  with  a 
frequency  of  a  few  MHz  is  strongly  absorbed  in 
the  lower  ionosphere,  resulting  in  an 
Instantaneous  increase  in  the  electron 
teMperature  T.  The  doalnant  energy  loss  process, 
excitation  of  rotational  Nj,  levels,  is  also  very 
fast  (i.e.  less  than  i  asec),  resulting  In 
Instantaneous  teaperature  relaxation,  when  the  HP 
power  is  reduced.  Transmission  of  lui  amplitude 
Modulated  HP  wave  results  In  oscillation  In  the 
electron  temperature  followed  by  corresponding 
oscillation  in  the  conductivity  tensor.  If  the 
HP  heating  occurs  In  ionospheric  regions  where 
aabient  d.c.  currents  flow,  an  a.c.  current  is 
generated  In  the  heated  volume  which  radiates  at 
the  modulation  frequency  CChang  et  al.  1981, 
Stubbs  et  al.  19811.  In  this  letter  we 
deaionstrate  that  the  presence  of  a  d.c.  current 
is  not  necessary  and  a  substantial  a.c.  current 
at  the  modulated  frequency  can  be  spontaneously 
generated  by  a  dynaeio  process.  The  efficiency  of 
the  process  is  Inversely  proportional  to  the  low 
frequency.  It  is  also  suggested  that  this 
process  rather  than  current  modulation  is  the 
source  of  the  PCS  pulsations  received  by  Stubbe 
and  Kopka  ( 1981 )  during  ionospheric  heating 
experiments. 


X  (7XB) 
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In  Eq.  (1)  n  is  the  electron  density,  u  is  the 
plasma  flow  velocity,  a  the  conductivity,  and  the 
aabient  magnetic  field  «ms  taken  as  Bq  =  e2Bo. 

The  source  term  on  the  r.h.s.  of  Eq.  (i) 
represents  a  dynamo  current  due  to  the  presence 
of  a  density  gradient  perpendicular  to  the 
temperature  gradient,  it  is  obvious  from  Eq.  (1) 
that  in  the  presence  of  a  stationary  density 
gradient  (Vn)  and  an  oscillatory  temperature  T 
the  spontaneously  generated  sbgnetic  field  will 
have  a  similar  tism  dependence  phase  shifted  by 
n/2;  therefore,  the  heated  region  will  act  as  a 
radiating  euntenna.  It  is  not  difficult  to  solve 
Eq.  (1)  including  the  diffusion  as  well  as  the 
Hall  and  convection  terms.  It  is,  however, 
obvious  that  for  a  low  B-plasma  (B  defined  as  the 
ratio  of  plasma  pressure  to  Magnetic  pressure) 
the  convective  term  can  be  neglected,  in 
addition,  since  this  letter  addresses  only  the 
feasibility  of  the  process  and  not  the  details  of 
the  radiation  pattern,  we  consider  the  simple 
case  Where  the  aabient  magnetic  field  Bq  i* 
parallel  to  the  density  gradient  Vnt  this 
corresponds  to  the  ejqbrismntal  configuration  in 
the  auroral  zones.  The  aziamthal  (Bs)  component 
of  Eq.  (1),  is  then  given  by 

=  Is 

%dtere  K  a  c^/t,ira  is  the  magnetic  diffusivlty  in 
the  aabient  plasab. 

In  order  to  provide  quantitative  estimates  of 
the  amplitude  Bg,  we  assume  cylindrical  symmetry 
about  the  z-direction  and  prescribe  the  source 
term  in  Eq.  (2)  as 

S(r,z,t)  ■  ^  (Vn)^(VT)^  = 


Generation  of  Oscillatory  Magnetic  Field  S  r  i  r  **  i 

=  —  expl - -]  expl - -j  [l-coe(u)t)J  (3a) 

Spontaneous  generation  of  oiagnetic  fields  froa  Ltp 

the  interaction  of  electrosiagnetic  (em)  waves  where 

with  plassbs  was  demonstrated  by  Stamper  et  al . 

(1971)  in  a  laser -plasma  interaction  experiment.  c  1  T 

A  staple  equation  was  derived  in  this  paper  for  ®o  "  e  iT  L~  ’  (3^) 

the  spontaneously  created  magnetic  field  B  in  the  n  T 

fora  of 


Iff  X  ((Vnle/nl'^  is  the  aabient  density  gradient 

Copyright  1989  by  the  American  Geophysical  Union.  length,  T/l/j  la  the  voliass  averaged  temperature 

gradient  in  the  HP  heated  plasma,  and  Lg 

paper  Nkaber  9L643S  represents  the  characteristic  dlawnsion  of  the 

O094-8278/8S/00SI.-843S803.0O  heated  region  in  the  vertical  direction. 
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The  cylindrically  sysnetric  version  of  Eq.  (2) 
with  the  source  tens  given  by  Eqs.  (3)  can  be 
solved  as  generalized  diffusion  equations  by 
using  a  coariblnatlon  of  Fourier -Bessel 
transfozsw.  Following  the  procedure  outlined  in 
the  appendix  (Eq.  AS)  we  find 

o 

r  r 

(i-cosT.)exp[-  g^^j,xp[- 

(D(t-t-))  (FfT-T-))** 


r=Mt ,  T • =ut •  and 


D<T-T'  ) 


P(T-T-  ) 


(T-T-  )+t 


(T-T*  )-H 


nie  aeplltude  of  the  spontaneously  generated 
■egnetic  field  Bs  will  be  controlled  by  the 
values  of  4X/uL>r*  and  .  These  parasMters 

correspond  to  the  ratio  of  the  radial  and  axial 
diffusion  tine  scales  and  to  the 

ELF/ULF  tisM  scale  4/u.  They  represent  the 
balancing  between  Mgnetic  field  dissipation  due 
to  the  diffusion  ten  in  Eq.  (2)  amd  the  nagnetic 
field  generation  due  to  the  source  ten. 

Equation  (4)  has  staple  solutions  in  the  limiting 
cases  of  4K/u(t^*.  Cz*)  »  i  and  4K/<j(i:/r*,  t«z*) 

«  i.  In  the  first  case  the  diffusion  tm 
through  the  modified  region  is  short  compared  to 
the  ELF/ULF  time  scale  (4/u)  and  limits  the 
amplitude  to 

3 

Bg(0,0,t)  =  Cmin(L^,L*)lCl-coe(ut)l  (6) 

Eq.  (6)  was  found  by  approximating  the  time 
integration  by  taking  the  most  significant 
contribution  from  the  integrand  at  t’=t  and 
■Kiltiplying  by  the  width  of  the  contributed  part 
tmin  (L<r^.Lz*)l  u/8K.  This  simply  states  the 
physically  obvious  conclusion  that  the  field 
amplitude  is  SoTdiff.  where  is  the  minimum 

axial  or  radial  diffusion  time  through  the 
andified  region. 

The  second  limit  corresponds  to  long 
diffusion  time  relative  to  the  EIT/VLF 

build  up  tism  4/u.  The  value  of  Bg  found  by 
performing  the  tism  integration  in  the  limit  that 
k«0  is 


B^(0,0,t)  =  —  tut-sin(ut)] 


ULF  Generation  in  the  Lower  Ionosphere 


From  Eq.  (7),  we  can  see  that  the 
spontaneously  generated  magnetic  field  consists 
of  two  parts.  A  steady  isagnetic  field  buildup 
increasing  linearly  with  tism  due  to  zero 
diffuslvlty,  and  an  oscillatory  piece  due  to 
temperature  smdulatlon.  The  amplitude  of  the 
oscillatory  smgnetic  field  is  Sq/Iu,  and  in 
practical  units,  will  be  given  by 


IB  (r.z.t)!  = 


-1  f  HzU  T  \  flOkml  fl 
[f  J|iOOO.xJ[Lj|- 


An  important  result  of  Eq.  (8)  is  the  scaling  of 
the  field  amplitude  as  1/f.  Stubbe  and  Kopka 
(1981)  detected  a  B  field  amplitude  of  10. 8Y  for 
10  min.  and  5.4y  for  5  min.  period  micro - 
pulsations  excited  by  heating,  which  is 
consistent  with  t)ie  l/f  dependence  of  Eq.  (8). 
Notice  that  if  we  tidee  T  «  10^,  L(|  »  Lp  «  10  km 
which  corresponds  to  expected  eiqmrinmntal 
values,  we  find  from  Eq.  (8)  Bg  «  5y  for  f  =-  .02 
Hz  (T  =  S  min. )  and  Bg  »  lOy  for  f  «  .01  Hz.  As 
mentioned  by  Stubbe  and  Kopka  (1981)  the  power 
expected  on  t)ie  basis  of  current  modulation  for 
the  PCS  frequencies  is  by  two  orders  of  magnitude 
smaller.  In  comparing  the  estimates  of  Eq.  (8) 
with  the  observations  on  t)>e  EtF-VLP  range  it  is 
clear  that  the  l/f  dependence  does  not  hold  for 
the  range  above  200  Hz  and  that  the  values 
expected  from  Eq.  (8)  are  much  lower  than  the 
ones  eiqiected  on  the  basis  of  current 
modulation.  For  example  for  f  =  2)(Hz,  Eq.  (8) 
predicts  lo'*y  at  the  source  vdilch  is  much 
smaller  than  the  value  of  io*’y  measured  on  the 
ground.  On  the  other  hand  the  values  detected 
near  and  below  10  Hz  are  consistent  with  the 
observed,  on  the  basis  of  the  above  we  expect 
that  the  transition  from  current  modulation  to 
dynamo  action  being  the  dominant  process  occurs 
in  the  range  between  10*200  Hz.  We  note  that 
inclusion  of  the  Hall  term  in  Eq.  (1)  does  not 
affect  the  total  magnetic  field  energy  produced 
by  dynamo  action;  the  Hall  term  simply 
redistributes  the  Bg  generated  flux  to  Bj.  and  B^ 
components.  For  the  E- region  and  Lp  “  10km  the 
long  diffusion  time  approximation  is  restricted 
to  frequenices  f  >  .OlHz.  Below  this  frequency 
field  saturation  described  by  Eq.  (6)  occurs. 

Power  Estimates 

From  Eq.  ( 7 )  we  find  that  the  power  radiated 
in  the  low  frequency  wave  is 


_  1  2  a  e 
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which  in  practical  units  becomes 


Pl^  *  .3 


10km  IHZ 
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In  this  letter  we  will  concentrate  on  the 
limit  of  long  diffusion  times,  which  for  the 
frequenices  under  consideration  is  applicable  in 
the  lower  ionosphere.  The  situation  described  by 
Eq.  (6)  as  well  as  intermediate  cases  will  be 
discussed  in  a  future  publication. 


Since  values  of  T  of  the  order  of  4-5  x  10’*K  can 
easily  be  obtained,  radiating  power  of  the  order 
of  several  tens  of  Watts  can  be  generated  in  the 
UtF  and  Icwer  EIF  range.  Such  temperatures 
require  powers  of  1-2  MH  CPerkins  and  Roble  1978; 
Chang  et  al.  19811.  Optimization  of  the  schesie 


Our  startins  ooint  for  the 
currents  that  couple  the  three 
Perkins  and  Goldman,  19813  the 
fluid  eouations 


derivation  of  the  nonlinear 
modes  is  LGure\rich,  1977; 
Marm  collisional  electron 
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Mhich  tosether  with  the  equation  of  state  for  the  pressure 
Oe  ~  ’^e  ^e'  form  a  closed  svs  .iem  of  equations  for  Og »  Vg  ► 
Tg’,  •  Av  is  the  mean  fraction  of  enersv  lost  per  collision, 
and  Tq  the  ambient  temperature.  In  (10)  we  have  neslected 
heat  conduction,  an  assumption  valid  for  the  E  but  not  for 
the  F  region.  Other  therms  neslected,  such  as  friction  and 
thermal  force,  do  not  chanse  the  final  results.  The 
nonlinear  currents  can  be  computed  by  considerins  small 
perturbations  about  a  nondriftins  equilibrium  (  Oq  ,  Tg)  so 


that  Vg 


V,  n-  =  no  + 


o^n  ‘ 


=  T g+T ,  assumins  plane  wave 


solutions  of  the  form  expC- i (wt-k . i ) 3  and  usins  the 
resonance  conditions 


‘^1  "  *^2  "  ^3 


(II) 


The  computation  is  tedious  but  straishtf onward  and  is  siven 
in  the  .appendix.  The  resultine  hiah-f  requency  nonline.ar 
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which  allows  propasation  of  Alfven  waves  as  Ions  as  v^/u^  << 

1. 

( 2 )  For  (j„  >  >  w  •  v_ 

3  Cl  3 

XX  ^ 

€3  =  0  (6a) 


J  u  (J_ 

ce  3 

which  IS  the  usual  electron  whistler. 


(3)  For  V2  >>  tJci*  ^3 


=  0 


3  <i» 

ce  3 

resardless  of  the  ratio.  This  is  similar  to  the 

whistler  dispersion  found  above  for  U3  >  Wci'  is  the 

helicon  mode  (Menyuk»  Ko#  and  Paoadoooulos*  manuscript  in 
preparation.  1984).  In  this  case,  contrary  to  case  (I) 
above,  the  current  is  carried  by  the  electrons  while  the 
ions  are  viscously  frozen.  In  this  sense  it  is  analosous  to 
the  well-known  helicon  wave  in  solid  state  physics  CAisrain. 
1961]  in  which  the  current  is  carried  by  the  free  electrons, 
while  the  ion  motion  is  frozen  by  the  lattice.  Referrine  to 
the  excitation  of  ELF  waves  below  lOQ  Hz  in  the  lower 
ionosphere  it  is  obvious  that  the  relevant  mode  is  the 
helicon  rather  than  the  Alfven  mode.  to  which  we  shall 
restrict  our  attention. 
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low-frequency  wave  is  then  CKrall  and  Tr  ivelp  iece ,  1973.1 


u 


3  XX'  ,  2 
?  '3  -  ''3^ 


(jJ 


xy 


u 


3  xy 
2  ®3 


(i) 


XX 


=  0 


(2  ) 


2  3  3 

c 

are  the  diasonal  and  the  off-dia«onal 
elements  of  the  cold  plasma  dielectric  tensor.  Assumins 


Xlf  V 

where  £3  and  £3 


Wee 

>  >  V  >  >  CJ3  we 

find 

exx 

£  =1 

2 

(1)  u., 

pe  3 

exy 

2 

(li 

pe 

(  3a) 

®3  ^ 

~T~  v~ 
CJ 

ce 

3 

(•>  (li¬ 
ce  3 

ixx 

2  , 
*^01  <"3 

>  1V3) 

i  xy_ 

.  2 

^“'pi  "ci 

(3b) 

*3 

"3["ci‘“'3 

+  ■•'31^ 

®3  “ 

"3 

"L"  ‘“3  *  ‘•'3’^. 

For  Upi  >  v.max  ((Jci’  *'3*  equations  (3)  reduce  l.o 

2 


XX  ,  exx  ixx 

'3  =  ‘  '3  *  '3 


^xy  _  ^exy  ^  ^ixy 
3  3  *3  ^ 


(V 


o  1 


U3+.V3 


.  -  (W_+1V»)^1 

Cl  3  3 


Cl/ 


OJ 


pe 


w  w_ 
ce  3 


1  - 


(CV3+1V3) 


(4a) 

(  4b) 


(j 


c  i 


-1 


From  (2)  and  (4)  and  noelectins  the  displacement  current  we 
find  the  followine  reeimes  of  weakly  attenuated  proper 


e leenmodes : 


compute  the  power  and  electric  field  threshold  conditions 
for  ELF  seneration  in  the  lower  ionosphere  by  two  HF  pumps 
under  steady  state  conditions  and  to  discuss  the  antenna 
requirements  for  a  proof  of  principle  experiment.  The  final 
section  summarizes  our  results,  stresses  their  limitations, 
and  outlines  potential  improvements  and  extensions. 

2.  NONLINEAR  INTERACTION  OF  TWO  HF  WAVE  PACKETS 
WITH  AN  ELF  WAVE  IN  THE  LOWER  IONOSPHERE 


We  assume  first  that  the  two  interactins  h i sh-f reauency 
waves  have  (jj.  Wj  >>  ‘**ce*  where  0)ce  electron 

cyclotron  frequency.  In  this  case  the  plasma  behaves 
isotropically,  so  that  the  dielectric  tensors  £1,3  of  the 
two  waves  are  diasonal.  i.e.,  €1,2  I*  fluid 

approximation  the  dispersion  relation  is  then 


O 


1.2 


u 


1.2 


2 

c 


e 


1.2 


(  la) 


€ 


1.2 


1 


2 

<t> 


pe 


(lb) 


where  (Jpe  is  the  plasma  frequency  and  vthe  electron-neutral 
collision  frequency.  For  the  low- frequency  wave  we  consider 
W3  <  Wpj.  where  Upi  is  the  ion  plasma  frequency,  so  that  we 
can  neslect  its  electric  field  parallel  to  the  ambient 
magnetic  field,  whose  direction  is  taken  alons  the  z  axis 
(i.e..  £3^  “  0)  .  The  dispersion  relation  for  the 
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Gurevicht  1978]  dominates  over  the  ponderomot i ve  force,  the 
downconvers ion  efficiency  is  increased  by  a  factor  v/u^ 
where  v  is.  the  electron  collision  frequency.  Furthermore 
the  power  threshold  required  for  complete  decay  is  reduced 
by  a  comparable  v/u)^  factor.  Since  the  valm;!;  of  v  are 
larser  for  the  lower  ionosphere  (i.e.,  E  resion),  it  is 
advantaseous  to  consider  the  interaction  in  this  resion. 

In  this  paper  we  examine  the  ELF  seneration  in  the 
ionospheric  E  resion  by  the  beatins  of  two  HF  heaters.  In 
the  next  section  we  discuss  the  ELF  plasma  eisenmodes  in  the 
E  resion  where  the  two  HF  pumps  will  couple  and  derive  the 
equations  describins  the  nonlinear  interaction  of  the  three 
wave  packets.  It  is  noted  that  for  altitudes  below  150  km 
the  lon-neutraL  collision  frequency  (V3)  is  larser  than  the 
ion  cyclotron  frequency  (Uci ) >  thereby  not  allowins  the 
existence  of  proper  Alfven  or  masnetosonic  eisenmodes.  It 
is  then  found  that  the  viscous  ion  dampins  (V3  >  Wci^ 
freezes  the  ion  motion  and  allows  the  whistler  mode  to 
continue  as  a  helicon  CAisrain,  1961]  a  risht-and  polarized 
mode,  to  frequencies  W3  <  V3  .  In  Section  3,  we  extend  the 
work  performed  on  threshold  conditions  and  power  transfer 
efficiency  in  col 1 i s 1 onless  three  wave  packet  interactions 
to  the  collisional  resime.  This  section  includes  many  novel 
and  rather  profound  physical  issues.  Section  4,  uses  the 
mathematical  and  numerical  conclusions  of  Section  3,  to 


nonlinear  mechanisms  for  ELF  seneration  which  are 
independent  of  the  ambient  current.  Paoadopoulos  et  al . 
C1982!)  considered  the  possibility  of  parametric  excitation 
of  an  Alfven  or  masnetosonic  wave  in  the  ionosphere  by 
beatins  two  high-frequency  waves  with  equal 
to  the  frequency  of  the  low-frequency  mode  (CJ3).  The 
process  envisioned  was  essentially  a  modified  stimulated 
forward  Brillouin  scatterine  CLiu>  19761  with  the  Alfven 
wave«  instead  of  an  ion  acoustic  wave*  actina  as  the 
low-frequency  signal .  In  the  language  of  parametric 
instabilities  the  pump  (u^)  and  the  idler  (oij)  are  HF  waves* 
while  the  signal  ((*>3)  was  an  Alfven  or  masnetosonic  mode. 
In  the  work  by  Papadopoulos  et  al.  119821  the  process  was 
considered  as  col 1 isionless  so  that  the  ponderomot ive  force 
was  the  dominant  nonlinear  force.  The  power  threshold 
necessary  for  complete  decay  can  be  found  by  the  inverse 
scattering  transform  (1ST)  method  ZZakharov  and  Manakov, 
1975;  Reiman*  1977;  Xaup  et  al.*  19791  The  aqalysis 
revealed  that  for  upper  F  region  ionospheric  parameters  the 
power  threshold  necessary  to  excite  Alfven  waves  was  very 
large.  In  addition  the  maximum  downconversion  efficiency 
was  limited  by  the  Manley-Rowe  relations  CSagdeev  and 
Galeev*  19691  to  W3/CJ1*  i.e.*  less  than  10  for 
downconvert  ins  5  MHz  to  100  Hz.  It  was  subsequently  noted 
CPapadopoulos  et  al.*  19831  that  if  the  thermal  nonlinearity 
caused  by  collisional  electron  heating  [Berger  et  al . *  1975; 


Doe  to  technical  problems  with  the  Tromso  heat  ins 


facility  it  has  not  been  possible  to  senerate  sisnals  in  the 
lO-Hz  to  2Q0--HZ  resime.  However*  sisnals  up  to  10  y  have 
been  senerated  for  Pc  5  pulsations  up  to  10  Hz.  These 
sisnals  are  much  larser  than  expected  from  current 
modulation*  and  their  .jhysical  orisin  has  been  a  mystery. 
In  a  recent  paper*  Papadopouios  and  Chans*  1985  noted  that 
spontaneously  senerated  masnetic  fields  with  time  variation 
similar  to  the  HF-  modulation  can  be  senerated  durins 
ionospheric  heatins  experiments  when  Vp  X  VTq  ^  0*  where  Vp 
is  the  ambient  ionospheric  density  sradient  and  VT^  the 
electron  temperature  sradient  in  the  heated  resion.  The 
senoration  of  these  fields  is  independent  of  the  ambient 
current*  and  its  physical  orisin  has  been  demonstrated  m 
laser-produced  plasmas  CStamoer  et  al.*  1971].  Papadopoulos 
and  Chans  11985]  indicate  that  for  frequencies  below  40  Hz 
the  spontaneous  field  seneration  dominates  over  the  current 
modulation*  thereby  explainins  the  larse  observed  sisnals. 

While  the  above  mechanisms  appear  to  be  rather  well 
understood*  they  are  relatively  inefficient  in 
downconvo rt i ns  HF  power  to  ELF.  In  addition  the  nonlinear 
demodulation  can  be  applied  only  in  resions  of  preexistins 
currents*  and  the  maximum  senerated  sisnal  amplitude  is 
strictly  constrained  by  the  local  current  values.  These 


limitations  provoked  the  search  for  other  more  efficient 


waves  amplitude  modulated  at  the  desired  low  frequency 


The 


basic  physics  of  the  technique^  known  as  nonlinear 
demodulation  or,  for  historical  reasons,  ionospheric 
detection,  lies  in  the  fact  that  the  HF  wave  produces 
modulated  electron  heatins,  which  results  in  modulation  of 
the  local  conductivity.  Natural  ionospheric  currents  which 
pass  throush  the  illuminated  resion  are  modulated  by  the 
conductivity  chanses,  producins  a  radiatins  dipole  current 
pattern  ZChana  et  al . ,  19811.  The  earliest  experimental 

results  on  nonlinear  demodulation  were  observed  durins  the 
Soviet  heatins  experiments  at  Gorki i  CGermantsev  et  ai . , 
1973;  Kapustin  et  al  .  ,  19773.  The  polar  electrojet  was 

modulated  in  the  VI. F  ranse  (1-7  kHz)  with  a  5.75-MHz  carrier 
and  15-MW  effective  radiative  power  (ERP).  Sisnals  were 
observed  only  between  2.5  and  4  kHz  with  strensth  in  the 
ranse  2-25  x  10"^  ^V/m.  More  reliable  results  for  ool;ir 
electrojet  modulation  were  reported  by  the  Max  Planck  sroup 
CShubbe  et  al.,  1982a,  b)  us  ins  the  Tromso  Norway 
ionospheric  heater.  The  ERP  in  this  case  varied  between  75 
and  125  MW,  and  the  HF  frequency  utilized  was  2.5-8  MHz. 
VLF  sisnals  with  amplitude  10“^  y  were  observed  with  maxima 
for  2  kHz.  Fin.ally  Farrnro  et  al.  C19821  reported  VLF 
sisnals  in  the  500-Hz  to  5-kHz  ranee  usins  the  Arecibo  HF 
heater  ,jt  .1.17  MHz  to  modulate  the  equatorial  electrojet. 
Asain  maximum  strensth  appeared  for  2-kHz  VLF  sisnals. 
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1.  introcxx:tion 

The  interaction  of  hish-power  radio  waves  with  the 
ionospheric  plasma  is  stronaly  nonlinear.  The  main 
nonl 1 near  1 1 1 es  arise  throush  the  ponderomot i ve  and  thermal 
forces.  As  a  result  of  the  nonlinear  interactions, 
electrostatic  (es)  and  electromaanet ic  (em)  waves  with 
frequencies  different  than  the  pump  frequency  can  he 
senerated.  A  comprehensive  review  of  the  subject  can  he 
found  in  a  recent  monosraph  LGurevicht  19783.  The  present 
paper  examines  the  feasibility  and  the  efficiency  of 
utilizins  the  nonl inear ities  of  the  ionospheric  plasma  to 
downconvert  HF  sisnals  to  ELF  in  a  controlled  fashion.  The 
subject  matter,  besides  its  intrinsic  scientific  merit,  has 
a  broad  spectrum  of  applications,  ransins  from  ionospheric 
and  masnetosphe r 1 c  probins  to  low-frequency  submarine 
comnun  1  cat  ions .  The  emphasis  in  this  work  will  be  on  the 
nonlinear  physics  aspects  of  the  interaction,  the  projected 
scalinss.  and  the  heater  requirements  for  a  proof  of 
principle  experiment. 

Ionospheric  heaters  have  been  used  successfully  for  the 
seneration  of  VLF  and  UL  F  sisnals  CStubhe  and  Kopka  1977; 
Stubbe  et  al  .  ,  1981;  Ferraro  et  al  .  ,  19823.  These 
experiments  induced  local  modulation  of  existins  ionospheric 
currents,  such  as  the  aiiroral  and  equatorial  electrojets,  hy 


illuminatins  the  cur rent-carryins  resion  with  strons  HF 


Generation  of  ELF  waves  by  stimulated  parametric  couplins 
of  two  HF  waves  in  the  lower  ionosphere  is  considered.  In 
this  resion  the  nonlinear  force  is  dominated  by  the  thermal 
rather  than  the  ponderomot i ve  nonlinearity.  It  is  shown 
that  this  results  in  lowerins  the  pump  threshold  for 
complete  decay  by  more  than  an  order  of  masnitude  while 
achievins  efficiencies  in  excess  of  those  expected  on  the 
basis  of  the  Manley-Rowe  relations.  The  lowei — frequency 
mode  excited  for  E  resion  altitudes  is  the  helicon  mode* 
which  continues  to  frequencies  below  the  ion  cyclotron 
frequency  because  ion-neutral  collisions  freeze  the  ion 
motion.  The  application  of  these  results  to  ELF  seneration 
in  the  lower  ionosphere*  includins  power  estimates  for  a 
proof  of  principle  experiment*  is  discussed. 
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roqulrM  Huciaising  the  electron  teeperaturee  (T* 
dependence)  %diile  ainislzing  the  heated  area 
(l/lfl  dependence).  Both  can  be  achieved  with 
high  antenna  gain  and  by  optlnlzlng  the  energy 
depoeltlon  altitude.  Detailed  nunwrical  atudlea 
are  In  progreaa  and  will  be  reported  elaeidiere. 

SuBMiy  and  Dlacuaalon 


idtere  Be(t,k.o)  la  the  Fourier -Beaael  coeponent 
of  tlte  initial  condition  Be(r,a,o)  at  t  s  o. 
Pinal  aolution  of  Zq.  (Al)  la  then  obtained  by 
the  Inverae  tranafon 

OD  CD 

®e(r,E,t)  a  J  dl  I  J  dK  Be(f.li.t)  Jo(lr)*^*^ 

O  -CD 


Ne  have  ahown  that  aaplltude  nodulated 
lonoapherlc  heating  can  generate  apontaneoua 
nagnetic  fielda  by  coupling  to  the  denalty 
gradlenta  of  the  lonoapherlc  plaaana.  The  tian 
dependence  of  the  apontaneoua  magnetic  fielda  la 
almilar  to  the  dependence  of  the  tenperature 
gradient  perpendicular  to  the  lonoapherlc 
tenperature  gradient,  ahlfted  by  ir/2  In  phaae. 
Therefore  a  periodic  heating  pulae  will  generate 
periodic  magnetic  pulaatlon.  nie  aaplltude  and 
power  of  t)ie  low  frequency  algnal  acalea  of  l/f 
favoring  t)te  lower  frequency  algnala.  It  la 
laportant  to  note  that  the  awchanlam  doea  not 
require  and  la  not  affected  by  the  preaence  of 
ambient  lonoapherlc  currenta.  nte  field 
generation  la  due  to  the  preaence  of  a  aource 
term  producing  fielda  with  V  x  B  ^  0  and  la  not  a 
plaama  or  parametric  inatablllty  auch  aa 
dlacuaaed  by  Kuo  and  Lee  (1983).  In  thle  paper 
we  daannetrated  the  baalc  phyalcal  mechanlam 
along  with  approximate  eatlmatea  of  the  ei^ected 
aaplltude  and  power.  It  la  eaay  to  ahow  that 
abort  acale  denalty  fluctuatlona  due  to 
ponderoantlve  force  or  denalty  eudulatlon  will 
not  affect  the  aero  order  field  generation.  The 
role  of  dlffualon  aa  well  ae  application  and 
mealing  of  the  proceaa  for  altuatlona  that  7n  la 
at  angle  to  Bo>  are  preaently  under  atudy  ^md 
will  be  reported  elaetdiere. 

appendix 

For  cylindrical  ayaawtry,  such  as  e7q)ected 
In  the  HF  heating  experlamnt,  Eq.  (2)  reduces  to 


Note  that  this  procedure  leaves  the  form  of 
source  term  S(r,E,t)  unspecified,  and  can  be 
extended  to  include  the  asimuthal  dependence  of 
the  problem  by  adding  a  e  transform 
J  de  eiN®, 

For  the  source  term  given  by  Eqs .  ( 3 )  and 
Be(r,a,o)xO,  sq.  (AS)  gives 
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there  t  =  ut,  t’  *  ut’  and  the  Integral  relations 
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I  dr  ^ 
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*  -iXz  -(*/L  _  -(L  V2)* 

da  e  e  s  /n  L^e  (A7 


were  used  to  obtain  the  final  result. 


r,z,t)  s  S(r,a,t) 
(Al) 


there  S(r,a,t)  la  an  arbitrary  time  dependent 
source  given  by  Eqs.  (3).  Introducing  a  Bessel 
transform  In  r,  and  a  Fourier  transform  In  a, 
such  as 


OD  00 

^(l,k,t)  »  1^  ^  dr  r  I  da  Bg(r,a,t)Jo(  lr)e 


-ika 


{A2) 


equation  (Al)  reduces  to 


I  It  *  Be(l,k,t)  =  S(l,k,t)  (A3) 


This  aquation  la  first  order  In  time  and  yields 
an  exact  solution 


B^d.k.t)  =  Bg(*,k,o)e 


(A4) 


♦  J  dfS(  l.k.f  )e‘^  *  )K(t-t') 
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currents  are 


<  1  - 


4wm  G)_ 
0  ^ 


)<k^ 
CJ^  3 


®3  ' 


(12) 


=  7^ - i  (1  -  — )(k  •  e®  .  E„)* 

2  4irm  (j,  (j  3  3 

e  1 


(tt  3 

Here  and  in  the  rest  of  the  caper,  v,  n,  T,  J,  and  E  are 
Fourier-transformed  quantities.  With  the  exception  of  the 
small  iv/a}i,2  terms,  (12)  is  the  standard  result  CWeiland 
and  yHhelmson,  19771,  The  low-frequency  nonlinear  current 
is  siven  by 


J^*-  =  -2_ 

“3  e 

*  k  (E  * 

3  1 

E*  )  (1 

4^  2 

J^)  (13) 

30,  3 

3  Anm 

a  "l“2  3 

For  v=0  this 

is  asain 

noth  ins 

more 

than 

the  standard 

low-frequency 

current 

caused 

by 

the 

collisionless 

ponderomot i ve 

force  CFapadopoulos 

et 

al .  , 

19821.  The 

coefficient  is  due  to  the  nonlinear  perturbation 
caused  by  Ohmic  heatine.  For  V/W3  >>  1  it  becomes  the 
dominant  term,  introducins  a  V/W3  enhancement  factor  in  the 
current  as  well  as  a  7r/2  phase  shift.  Several  physical 
comnents  on  the  derivation  of  (12)  and  (13)  can  be  found  in 
the  appendix. 


Usins  the  currents  J^*-  as  a  source  in  the  wave  equation 
for  each  of  the  three  waves  we  find 
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In  the  absence  of  nonlinear  currents  (J™*-  =  0)  the  vanishine 
of  the  determinants*  IDjI  =  0.  represents  the  linear 
dispersion  relations  of  the  individual  Maves*  that  were 
siven  by  (1)  and  <2).  Takins  without  loss  of  senerality  the 
polarizations  of  as 


^1.2  -^1.2^ 

and  usins  the  isotropic  property  of  the  tensors  61,2 
“i,2  >>  “ce,  we  find 


°1^1 


V2 


'*”^1  ,NL 
c 

‘■“2  ,NL 

~~r  ^2 

c 


with  Di 


in  by  <1)  and  2  (12>.  For  the 


polarization  siven  by  (15)  and  usins  the  fact  that  £3^  0 
for  W3  <<  Wp i »  the  equations  for  the  low-f reouency  wave  are 
siven  by 


F  +  -  - 

c 


xy 

3  3y 


6^'''  E 
3  3x 


XX  ,  2 

‘3  -  '‘3 


XX  xy 

where  and  are  siven  by  Eqs.  (4). 

from  (17)»  we  find 


‘  ‘”"3  ,NL  ,  ,  , 

- —  ‘'3X 

c 


'”*^3  IML 

J3^  (17b) 

c 


Eliminatins  £, 


'^3x  ^  ^3y 


where 


°3  = 


2  ®3 


3  XX  ,2 
1  ‘3  -  "3^ 


"3 

2 

c 

XX 

®3  - 

Xal 

to  Alfven 

waves 

or 

€3^^. 

will 

be 

9  are 

interested 

in 

couplina  to  heItcons>  so  that 

163^^1  /  €3^^  -  (nii/mg)  (V/Ugg)  >>  1. 

Therefore.  D3  becomes 


°3  =  - 


2  ,  ,  2 

*^3  xy  *^3 

k  k  _ -  I  €  ^  I  k  k  4-  — 

^3  2  '®3  '  3'^3z  2 

c  J  I  c 


with  €3^  aiven  by  (7b).  The  case  of  Alfven  waves  can  be 
examined  in  an  analosous  manner.  Assumins  without  loss  of 
senerality  that  k^y  =  0  we  find  that  k3  *63  •  E  ~  k^jf  €3^^ 
Ey,  so  that  J3  =  J^y  y.  From  (12),  (13),  (16),  (18).  and 
(20)  we  find  the  couoled  system  of  eauations. 


^  ^  0  (1 


LJ^  ,  !i!3 

*^2  "2 


/9  (1  -  — 


•r 

—  )  _L  _2. 

"1  *^1  *^3 


(21a) 


(2lb) 


Mhere 


3x  pe 


(21c) 


and  we  have  set  Eay  =  £3.  In  the  collisionless  case  the 
hi«h-f reouency  waves  produce  a  low-f reouency  ponderomot 1 ve 
force.  -eE^,  on  the  electrons,  which  results  in  a  nonlinear 
current  that  drives  E3.  In  the  collisional  case.  J5"-  - 
predominantly  doe  to  a  temperature  perturbation  when  the 
electrons  are  Ohmically  heated.  In  either  case,  the  beatine 
of  the  h  i8h-f  reouency  velocities  i ,  v^2  with  the  low- 
froouency  density  perturbation  nj*-  provides  the  primary 
couplins  mechanism  between  £3  and  the  h ish-f reouency  fields 
£3.  £3,  It  IS  worth  notins  that  for  interaction  with  U3  << 
Wj.  (J3 .  1^3!  <<  Ikil,lk2l>  and  by  choosins  the  HF  electric 
field  polarization  in  the  y  direciton.  we  end  up  with  all  of 
the  nonlinear  currents  and  the  low-f reouency  electric  field 
in  the  y  direction. 


Eouations  (21)  describe  the  interaction  of  infinite, 
homoseneous.  plane  waves.  In  practice,  however,  we  deal 
with  two  HF  wave  packets  interactins  with  an  ELF  wave  packet 
within  a  spatially  localized  resion  controlled  by  the 
antenna  sain.  In  order  to  determine  the  efficiency  of  the 
interaction  the  localized  and  1  nhomoseneous  nature  of  ttie 


Maveoacket  interaction  should  be  included.  For 


a  weakly 


nonlinear  interaction  between  localized  wave  packets*  we 
follow  the  •  standard  procedure  and  expand  Dj  in  a  Taylor 
series  about  the  phase  matchins  point  where  the  resonance 
conditions  (11)  are  met:  that  is.  we  let 


D  .  f 

j 


iao. 

•j 

d(ii 


It  *  '■j 


lao 

j 

ak 

j 


V  + 


operate  on 


E  (r,t) 

J 


11  a  .  (  r.  t )  e 

.7  J 


(  23a) 


where 

u  =  (J  /  (ao  /aw  .)^  (23b) 

J  3  3  3 


and  a  j 

is  the 

slowly 

varyina  envelope 

of 

the 

wav0 

packet . 

The 

1 1  near 

■lampina  rates  Fj 

are 

aiven  by 

the 

imaainary  parts  of  Oj  while  the  partial  derivatives  of  Dj 
are  evaluated  for  a  particular  mode  (bij.kj).  For  positive 
eneray  waves  as  is  the  case  here,  dDj/duij  >  0.  Usina  the 
fact  that.  to  lowest  order.  the  waves  obey  the  linear 
dispersion  relations  Dj  =  0.  we  obtain  the  equations  for  the 
collisional  three-wave  couplina  in  standard  foim 


(It 

^1 

•  V 

i  V  1 

]  ^2‘*3 

(It 

“2 

’  V 

^2) 

^2 

w  1 

=  -K  j 

r,  iv  1  » 

1  -  -  a,  a„ 

(j^J  13 

(It  ^ 

^^3 

.  V 

+ 

^3] 

^3 

1 

1! 

f,  ^  i4v  1  » 

1‘  ‘  30.3  )  ^1=3 

where 

(24a) 

(24b) 

(  24c) 
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=  - 


dai  .  die 
J  J 


(25a) 


are  the  sroup  velocities  of  the  waves  and 


K  -1/9/ 


ao^  ao^  ao  3 

d(jj^  d<J^  awg 


(  25b) 


is  the  syimietric  couplina  coefficient  for  the  col  1 1 s ionless 
interaction.  By  settina  d/dt  =  0.  one  finds  the  steady 
state  tMo-d imensional  version 


ll?  "  -1  I?  "  ^1  "l  =  %  1  -2-3 

(I?  ""2  §3?  ^2)  ""2  "  ~  ^1^3*^ 

[3?  ^3  §3r  ^  ^3!  ^3  =  ■"  ^1^2  ' 


^2^3 


3  T, 

[3?  ■"  ""2  33?  ^ 


(  26a) 


(  26b) 


(  26c) 


where  the  interaction  is  assumed  to  evolve  in  z  lujz  >  0) 
The  new  variables  are  expressed  as  follows: 


a-au  u=u  /  u  r=T/u  ( 27a ) 

J  J  JZ  j  jx  jz  j  j  jz 


K  =  K  /  (u,  u_  u,  ) 
Iz  ?.z  3z 


(  27b) 


Equations  (26)  with  the  definitions  (23)  and  (25) »  the  value 
of  0  aiven  by  (22).  and  Dj.  Dj’  O3  by  (1)  and  (20)  form  the 
basic  equations  of  our  problem.  In  the  next  section  we  will 
use  them  to  determine  the  transfer  efficiency  durina  the 
steady  state  interaction  of  the  three  wave  packets. 
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THRESHOLD  CONDITIONS  AND  TRANSFER  EFFICIENCY  FOR  THREE 


WAVEPACKET  INTERACTIONS  IN  DISSIPATIVE  SYSTEMS 

Before  attemotins  to  solve  Eqs.  (26)  for  the  situation 
under  study*  it  is  worthwhile  to  review  briefly  the  previous 
work  on  the  subject*  which*  due  to  its  rather  mathematical 
nature*  misht  not  have  reached  the  ionospheric  community. 

In  the  collisionless  case  (i.e.  v=0)*  the  couplins 

coefficients  are  symmetric  and  the  set  of  Eas.  (24)  has  an 
infinite  set  of  invariants.  This  allows  for  the  use  of 
inverse  scatterins  transform  techniques  (1ST)  (AbLowitz  et 
al.*  19B1).  The  one  dimensional  version  of  Eqs.  (24) 

(Zakharov  and  Manakov*  1975;  Reiman*  1977;  Kaup  et  al.* 

1979)  has  been  extensively  studied  in  the  literature  usins 
1ST  and  the  results  were  confirmed  by  numerical  solutions. 
As  noted  before*  the  two  dimensional  steady  state  Eqs.  (26) 
are  equivalent  to  the  transformed  one  dimensional  soace-time 
equations.  Therefore  the  results  of  the  one  dimension.al 
space  time  problem  can  be  directly  applied  to  our  practic.aL 
problem.  An  important  result  of  the  one  dimension.al 
space-time  collisionless  analysis  is  the  fact  that  the 
threshold  value  required  for  pump  depletion  is  equivalent  to 
the  pump  amplitude  reouired  to  achieve  an  absolute 
instability  within  its  width  L  (Reiman*  19/7). 

For  collisionless  throe  wave  interactions*  the 


eisenvalue  problem  Mhich  determines  the  conditions  for 


absolute  instability  in  a  rectansular  pumo  is  well  known 
(i.e.  the  backward  wave  oscillator)  (Bobroff  and  Haus  1967). 

^  Takina  the  pump  i nhomoaene i ty  to  be  alons  x»  the  instability 

is  absolute  if 

u_  u„  <0  ,  (28a) 

2x  3x 

w 

L  >  ^  L  .  (28b) 

2  c 

where  =  Iu2x^3xl^  ^  l')'ol  is  tbe  critical  width  and  I  >0  *  “ 
IK^ail  IS  the  uniform  medium  arowth  rate.  Eq.(28a)  requires 
the  decay  waves  to  be  oppositely  trav  .ina*  while  Eo.(28b)  ’ 
defines  the  pump  threshold  for  the  absolute  instability  to 
occur  within  a  width  L.  We  demonstrate  below  that  for  the 
stronaly  coiiisional  case  ( v/Ci;3  >>  1)  similar  considerations 
apply.  Note*  however*  that  since  the  couplina  coefficients 
are  by  a  factor  v/cjj  laraer*  the  required  pump  levels  are 
smaller  by  the  same  factor.  In  addition  the  transfer 

efficiency  increases  by  the  same  factor.  Notice  that  due  to 
the  dissipative  nature  of  the  system  1ST  techniques  are  not 
applicable  so  that  one  has  to  resort  to  numerical 

techniques . 


The  linear  equations  for  the  eisenvalue  problem  that 
determine  absolute  instability  when  one  considers  t)ie 
initial  stase  of  the  interaction  durins  which  the  pump  is 
larae*  a^  >>  33,3*  are  siven  by  the  one  dimensional  version 
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of  Eqs.  (24b)  and  (2^c).  Here  only  the  evolution  of  the 
decay  waves  is  important*  so  that  we  have 


lat  *  ax  *  ^2] 

[st  *  '^axh*  ^al  ^3 

Note  that  yo(.x)  contains 
def ined  by 


=  -y  ( X ) 
o 


-  -y  ( X ) 
o 


the  pump 


inhomoeene 1 ty 


(29a) 

{29b) 

and  A  is 


A  =  1  + 


i  4  V 


(30) 


For  the  threshold  calculation  hare*'  we  can  ianore  the  small 
collisional  correction  on  the  RHS  of  Eo .  (29a).  The 
linearly  couoled  system  is  solved  by  takina  the  conjuaate  of 
Eq .  (29a)  and  Laplace  transformina  in  time  as  e^^.  We  then 
make  the  substitution 


(a^  ,  a^ )  =  exp 


to  obtain 


1_ 

2 


fp+r.,  p+r. 


u 


2x 


u 


3x 


(A^.  A3) 


(31 ) 


where 


y  ( X) 

o 


u. 


A. 


(  32a ) 


2X 


[h  ~  ^3 


y  ( x) 

o 


u 


3x 


A2  A 


(32b) 


T) 

Equations  (  32 ) > 
conditions*  form 


2 


D+r, 


p-r3i 


u. 


(  33  ) 


2x  “3xJ 

when  subject  to  the  appropriate  boundary 


the  eisenvalue  problem  for  the  arowth  rate 
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p  which  may  in  seneral  be  complex 


t  =  L  /  lAI^  .  (35) 
c  c 

Equation  (34)  readily  admits  solutions  in  terms  of 
tr iBonomotr ic  functions  to  which  we  apply  the  boundary 
conditions  that  be  zero  at  -L/Z  and  A3  be  zero  at  L/2  for 
Ujx  ^  ^  This  leads  to  the  dispersion  relation 

cot  (o'*  L/1  )  +  1)1  =  0  ,  (36) 
c  c 

for  the  normalized  srowth  rate  i)ic.  Me  have  evaluated  (36) 
numerically  in  the  limit  V/W3  >>  1  which  sets  B  ~  v/2  and 
lAI  ~  v/U)^.  Complex  eisenvalues  are  found#  which  indicates 
that  absolute  instability  indeed  occurs  in  the  collisional 


res ime . 


The  real  and  imasinary  parts  of  ijl^  for  the  first 


four  sroMins  inodes  are  shown  in  Fisures  la  and  lb  as  a 


function  of  the  normalized  width  L/Iq.  For  comoarison,  the 
results  for  the  collisionless  resime  (v  =  0:  0  =  0,  A  =  1, 
Ic  Lc)  are  also  plotted:  in  this  case*  the  growth  rates 
are  purely  real.  From  Figure  la.  one  finds  the  threshold 
condition  in  the  collisional  regime  to  be 

/-  >  1.8  1  (37) 
c 

For  a  given  width  L.  it  is  clear  that  the  pump  value  a^ 
required  to  achieve  absolute  instability  for  the  collisional 
interaction,  is  smaller  by  the  factor  Iq/Lc  =  IAI~^.  Thus 
a  reduction  in  the  required  pump  power  Pi  *  I  ail*  by  the 
factor  W3/V  is  achieved;  in  the  lower  ionosphere,  this  can 
be  more  than  two  orders  of  magnitude.  The  specific  pump 
requirements  for  a  steady  state  interaction  relevant  to  ELF 
generation  will  be  discussed  in  a  later  section. 

As  noted  earlier,  a  key  result  due  to  Reiman  (1977)  for 
col 1 i sionless  interactions  is  that  if  the  threshold 
conditions  (28)  are  satisfied,  the  linear  stage  of  the 
absolute  instability  will  be  followed  by  pump  depletion  in 
the  nonlinear  stage.  It  has  been  shown  that  an  absolute 
instability  can  also  be  established  in  the  collisional 
regime.  Furthermore,  the  pump  threshold  is  much  lower  than 
in  the  col 1 i s 1 onless  case.  This  threshold  reduction  is  of 
great  practical  interest,  especially  for  applications  to  ELF 
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I 


senerat ion 


HoMever.  an  equally  imoortant  issue  is  whether 


the  absolute  instability  would  lead  to  pump  depletion  and 
if  so.  what  values  of  downconver  s  ion  efficiency  can  be 
obtained.  Me  demonstrate  below  that  nonlinear  saturation  of 
absolute  instabilities  by  pump  depletion  indeed  occurs  for 
collisional  interactions  and  that  downconversion 
efficiencies  hisher  than  those  predicted  by  the  Manley-Rowe 
relations  can  be  reached. 

The  feature  that  distinauishes  the  collisional  from 
the  collisionless  interaction  is  the  presence  of  dissipative 
terms  in  the  couplins  coefficients  so  that  the  system's 
action  and  eneray  are  no  lonaer  time  invariant.  To  see 
this,  we  obtain  from  (24)  the  followina  relations  for  the 
one  dimensional  space-time  evolution 


d 

Tt  ^12 

2  2  4v  *  1 

r,  la  J^+r„la-f  —  IK'IReCa,  a_  a.,) 

L  L  ^  c.  ^  ^  J  » 

(38a) 

a 

at 

OP 

.-J., 

-op 

rilajl^r3la3l^-  ^  IKI«.<a’  a^  33)] 

( 38b) 

a 

at  T 

*  5  vlKIRe<a'  a^  a3)] 

( 38c) 

where 

dx( la^ 1^ 

^  '=2'^  1-  °13'_J  ‘"'['^1'^  *  '»3'^1 

(39a) 

and 


21 


T 


3 

E 

j=l 


J  J 


dx  (dial 
J  J 


(39b) 


2 


Notice  that  the  value  of  /  dxlal^  is  nothins  more  than  the 
wave  action  <i.e..  number  of  quanta).  In  the  col 1 i s i onless 
limit*  the  sums  c^j  and  0^3  of  the  actions  and  the  total 
enersy  7  are  time  invariant.  Equations  (38)  with  zero  on 
the  riaht  hand  side  are  commonly  referred  to  as  the 
Manley-Rowe  relations.  A  direct  consequence  of  these 
relations  is  the  limitation  in  the  downconversion  efficiency 
from  the  pump  (Ui)  to  the  sienal  (Uy)  to  a  value  smaller  or 
equal  to  (jj^/oti.  Extensions  of  these  results  to  the  case  of 
weakly  dissioative  systems  with  ponderomot ive  nonlinearity 
(i.e.»  Fj/edj  <<  1)  rely  on  the  approximate  conservation  of 
Ci2*  Ci3*  and  7  for  them  to  be  solvable  by  1ST.  The  1ST 
results  have  been  compared  with  numerical  solutions  of  Eos. 
(24)  with  V  =  0  and  found  to  asree.  They  sive  the  same 
(jy/Ui  downconversion  efficiency  (Liu»  1976).  In  order  to 
examine  the  effect  of  the  thermal  nonlinearity  we  solve 
numerically  the  system  of  (24)  with  v  ^  0  but  Fj  =  0. 

For  the  sake  of  simplicity,  we  limit  ourselves  to  a 
model  problem  in  one  dimensional  soace-timo  (x.t).  As  a 
comparison,  the  results  for  the  coll isionless  interaction 
with  symmetric  couolins  coefficients  are  also  presented. 

s  are  chosen  accordins  to  (11): 
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The  frequencies 


of  the  wave 


for  numerical  convenience,  we  pick  in  arbitrary  units:  ~ 

10.1.  (t)2  ~  10,0.  V  =  1.2  and  Wj  =  0.1.  The  initial  wave 
envelopes  are  rectansular  and  the  decay  waves  have  small 
amplitudes:  specifically,  we  let  (  a^  I  /  I  a  ^  I  =0 . 01 .  Ia3l=:0.0. 

Because  we  have  in  mind  an  absolute  instability,  we  also  let 
the  decay  waves  to  have  opposite  velocities  "  ~^2x’ 

thus  satisfyina  (28a).  The  problem  is  solved  in  the  pump 
reference  frame  so  that  i-iix  '  pump  threshold 

condition.  (28b)  or  (37).  is  met  by  settina  la^l  =  1.0  in 
the  CO 1 1 1 s ion  less  case:  whereas,  in  the  collisional  case. 
lail  =  IAI~’*  =(4w/3  013)”**  =  0.25.  With  these  values.  the 
respective  normalized  widths  are  the  same.  i.e..  L/Lq  L/i^ 
=  3.2.  so  that  from  Fiaure  la.  one  expects  a  arowina  mode  in 
both  interactions. 

Usina  the  above  initial  conditions,  the  coupled  system 
of  eouations  (24)  is  numerically  intearated  forward  in  time. 
Fiaure  2  shows  the  t ime-asymptot ic  behavior  at  t  =  25.0.  We 
observe  tfiat  pump  depletion  occurs  in  both  interactions. 
In  the  collisonless  case.  the  decay  waves  emerae  as 
symmetric  pulses.  They  are  identified  as  solitons  as 
predicted  by  the  ISf  solutions  {Kaup  et  aJ  .  .  1979).  Ttiey 

are  consistent  with  the  Manley  Rowe  relations  which  predict 
an  asymptotic  state  with  Ia2l/la3|  ••  1  for  pump  depletion. 
On  the  other  hand.  the  collisional  interaction  exhibits 
nonsymmetric  pulses  with  33  >>  33.  which  suaaests  a  rather 


By  choosina  the  pump  and  sideband  field  polarizations  as 
Ei,2  =  El, 2  '**®  see  from  the  left  hand  side  of  (A, 2)  ii  j 
(A. 3)  the  associated  linear  density  and  temperature 
perturbations  vanish  (in  the  (x,z)  seometry). 
1 .  e. , 


since  the  h i ah- f requency  waves  are  elect romaanet ic .  For  the 
low-  frequency  wave,  the  maanetized,  cold,  linear  response 


as  aiven  by  the  left  hand  sides  of  (Al)  -  (A. 3)  is 


At  the  next  order,  we  evaluate  the  thermal  and 
nonlinear  terms  on  the  riaht  hand  sides  of  (A1)-(A.3),  usina 
the  lowest-order  solutions,  <A.4)-(A.8).  It  can  easily  be 
seen  that  there  is  no  nonlinear  modification  to  the  density 
at  all  three  frequencies  (i.e,  riaht  hand  sides  of  (A?) 
vanishes) , 
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APPEMOIX 


Tha  plane  wave  perturbations  in  a  resonant  three  wave 
interaction  satifsy  the  fluid  eauations  (8)-(10)  as  follows 
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The  subscripts  ut,  k  for  the  square  brackets  indicate  that 
the  appropriate  convolution  product  of  the  other  two  waves 
IS  taken  to  meet  the  matchins  conditions  of  (11).  Equations 
(A.1)-(A.3)  can  be  solved  perturbat i vely . 


Our  main  interest  concerns  the  decay  of  a  hi»h- 
frequency  e lect romasnet ic  pumo-wave  (Wi.  kj)  into  a  hiah- 
frequency  sideband  kj)  and  a  low  frequency  mode  (0*3  . 

kj)  in  the  collisional  resime.  If  we  assume  that  the  pump 
wave  and  its  sideband  have  frequencies  (113,2  ^ce‘ 

electrons  are  unmaanetized>  and  the  left  hand  side  of  (A) 
yields*  to  lowest  order,  their  driven  velocities  in  the 
linear*  cold  limit 
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electron  heatins  Mill  be  required  to  support  proper 
eMper imental  effort.  MeanMhile  the  paper  should  be  used  as 
a  suide  to. the  physics  expected  durins  the  interaction  and 
the  approximate  desisn  of  ionospheric  heaters  required  for  a 
proof  of  principle  experiment. 

Before  closins  we  should  note  that  S.  Ganauly  and  M. 
Gordon  (1984)  reported  preliminary  evidence  for  beat  Mave 
aeneration  for  frequencies  between  10  and  40  Hz>  usina  the 
Arecibo  HF  facility.  Their  results  were  consistent  with 
power  threshold  scalina  increasina  linearly  with  fj.  The 
preliminary  nature  of  ;;he  experimental  results  and  the  fact 
that  the  interaction  reaion  was  in  the  F  rather  than  the  E 
reeion*  where  heat  conditions*  nealected  in  our  analysis* 
dominate  the  eneray  transport*  does  not  allow  for  a 
quantitative  comparison  of  our  theory  with  the  experiment. 


5.  SUMMARY  AND  CONCLUSIONS 


Me  have  examined  stimulated  excitation  of  ELF  waves  in 
the  lower  ionosphere  by  the  use  of  two  HF  pumps.  Previous 
work  CPapadopouIos  et  a.1  .  ,  1982*  19831  addressed  the 
excitation  of  Alfven  and  masnetosonic  modes  by  similar 
techniques.  However,  since  for  ionospheric  conditions  Vj  ^ 
Uci  1-50  km.  excitation  of  these  modes  is  possible  only 
above  that  heisht.  The  thermal  ponderomot ive  force  drivins 
the  interaction  is  proportional  to  the  electron-neutral 
collision  frequency  v  which  is  by  more  than  two  orders  of 
masnitude  smaller  above  150  km  than  near  100  km.  The 
discovery  that  the  helicon  mode  CMenyuk  and  Papadopoulos . 
1984)  is  a  proper  eieenmode  of  the  system  even  for 
frequencies  with  <  V3  allowed  us  to  reduce  the  threshold 
for  excitation  of  ELF  waves  below  100  Hz  by  more  than  two 
orders  of  masnitude.  by  couplins  to  the  helicon  mode  in  the 
vicinity  of  100  km  altitude.  In  addition  to  focusins  on  the 
helicon  mode,  a  more  detailed  description  was  presented  of 
the  parametric  decay  processes  in  dissipative  media 
discussed  previously  in  a  letter  CPapadopouIos  et  al  .  . 
19831.  As  noted,  a  more  comprehensive  analysis  includins 
the  effects  of  oblique  propasation  of  the  pump  from  t)ie 
sround  to  the  ionosphere,  the  1 nhomoseneous  structure  of  the 
pump  in  the  interaction  resion.  the  i nhomoseneous  structure 
of  the  ionosphere,  and  a  self-consistent  description  of  the 


from  (50),  by  takins  k2x^^3x  "*  *5,  Tg  ••  10^  *K,  and  Ui/Upg 
•»  2.  The  value  of  "*  2  imolies  a  30*  incidence  angle 
of  the  transmitter,  which  is  sufficient  to  achieve  long 
propagation  paths  (Ginzburg,  1970).  Using  these  values  (50) 
becomes 


"thr  =  (^1  (^1  (tS^) 


(  51 ) 


o 

Table  1  shows  the  power  required  to  excite  100-Hz  and  50-Hz 
waves  as  a  function  of  height  as  computed  by  (51).  The  last 
column  shows  the  required  pump  frequency  versus  height  based 
on  30*  incidence.  Daytime  conditions  were  assumed  in 
compiling  Table  1  (Gurevich,  1978).  Notice  that  the  optimum 
interaction  height  is  100-km,  and  relatively  modest  power  is 
required.  Me  would  like  bo  stress  that  the  values  in  Table 
1  are  indicative  rather  than  exact.  We  feel  that  the 
assumption  of  Tg  10^  *K  is  very  conservative  and  Tg  3-4 
X  10^  *K  IS  more  realistic.  This  will  reduce  the  power 
requirement  by  factors  of  3-4.  On  the  other  hand  a 
variety  of  propagation  losses  have  been  neglected.  These 
will  be  considered  in  the  future.  It  is  interesting  to  note 
that  PthR  ”  (I/Hq)  (1/N).  Therefore  further  reduction  in 
^THR  can  be  achieved  if  we  could  increase  the  local 
ionization  by  long  pulse  radiation  previous  to  the  beat 
excitation. 


32 


we  note  that  the  ion  neutral  collision  frequency  V3  >  Ud  •" 
220  s~^  for  altitudes  below  130  km*  which  is  the  resion  of 
interest  in  the  present  work.  Therefore  to  excite 
loM-f requency  waves  below  100  Hz,  we  must  couple  nonl inear ly 
to  the  newly  discovered  helicon  branch  siven  by  Eos.  (7). 
As  discussed  in  Section  2*  this  branch  occupies  the  ranee  Uj 
<  Vj  and  extends  to  very  low  frequencies.  It  is  important 
to  notice  that  contrary  to  the  conventional  situation*  where 
the  currents  associated  with  the  low-frequency  waves  are 
carried  by  the  ions*  in  the  helicon  branch  the  currents  are 
carried  by  the  electrons.  The  ion  dynamics  is  viscously 
frozen  due  to  the  hiah  col  1  i s ional i ty  (V3  >  Wcj  )  .  Under 
these  circumstances  the  Alfven  and  maanetosonic  modes  are 
not  proper  eisenmodes. 

Fiaure  5  is  a  schematic  of  the  aeometry  in  which  two  HF 
waves  (cJi,  <013*  k2)  interact  nonlinearly  in  the 
ionosphere  to  drive  a  helicon  with  (J3  <  V3.  The  HF  pump  is 
incident  at  a  small  anale  to  0q  with  a  beam  width  equal  to 
or  smaller  than  the  local  density  scalelenath  Lp.  For  * 
U2  >>  W3  the  resonance  conditions  (11)  yield  a  k3  that 
propaaates  almost  directly  downward.  Such  an  interaction 
conf iaurat ion  can  be  achieved  by  arazins  incidence  of  the  HF 
beams  at  the  desired  ionospheric  heiaht  (Fiaure  6). 

Approximate  values  of  the  required  power  can  be  found 
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10^"k 


100  Hz 


MM  . 


The  threshold  power  scales  linearly  with  frequency «  so  that 
it  is  smaller  for  the  lowest  frequencies  for  which  the 
analysis  is  applicable.  The  optimum  interaction  resion 
heisht  can  be  found  by  maximizins  the  product  This  is 
due  to  the  weak  dependence  of  the  terms  (Ui/b/pg)  and 

exact  plasma  parameters.  Me  will  further 
elaborate  on  these  later  on. 


Me  proceed  next  to  use  the  above  equations  in  the 
desian  of  a  proof  of  principle  ionospheric  heatina 
experiment.  For  convenience  our  estimates  refer  to  the  case 
where  the  masnetic  field  is  perpendicular  to  the  density 
eradient  (Fisure  5).  This  misht  not  be  the  optimal  case. 
Determination  of  the  optimal  seometry  requires 
seneral izat ion  of  our  calculation  to  include  the  effects  of 
pump  propasation  between  the  around  and  the  interaction 
reaion*  and  the  i nhomoaeneous  structure  of  the  pump  near  the 
reflection  surface.  These  issues  are  currently  under  study 
and  will  be  reported  elsewhere.  Before  comparina  the 
threshold  conditions  for  the  nonlinear  interaction  we 
examine  the  type  of  low-frequency  waves  of  interest  for  our 
particular  application.  The  dispersion  relation  for  the 
low— f requency  ionospheric  modes  was  discussed  in  Section  2. 
Referrina  to  typical  ionospheric  conditions  (Gurevich,  197tt) 


(46) 


and  substitutins  for  i9  explicitly*  Me  obtain 

.r-  ,  1  ce  2x  L  1 

IE  (z=0)l  =  — —  5 —  -r XT-  TTV  I  .  . 

^  ®  (j2  '‘ax  ll 

pe  c 


The  requirement  for  absolute  instability  leadins  to  complete 
decay  is  siven  by  <42b).  Therefore  the  threshold  electric 
field  found  from  (46)  by  takine  L/Ic  2  is 


IE^(2-0) 


=  4  X  io2  (i?i2!|  ^  ‘  (47: 

'  *-  '  11?^)  “ax  lAi'  "■ 


For  ionospheric  conditions  Wce  =  7.5  X  10*  and  v  =  2.5 
X  lO”*-*-  N  Te  based  on  ion  neutral  collisions;  N  is  in 
cm“^  Usina  these  values  (47)  can  be  written 
in  practical  units  as 


(Ei(z=0)Ith«  =  22 


in5  -3 
10  cm 


2x  10  km 
...  L 


l“oeJ  rax 


N  T 


^3  I  v 
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The  power  required  to  achieve  this  can  be  found  from  (47)  by 
assumina  that  the  HF  illluminated  reaion  is  a  square  with 
side  L.  Then 

"tHR  “  ?  ‘o  ■= 

Notice  that  the  result ina  equation  is  independent  of  L. 
From  (48)  and  (49)  we  find  that  in  practical  units 


case  apply  to  the  steady  state  interaction  as  well.  Rather 
than  enersy*  here  one  looks  at  the  downconversion  of  HF 
power  into  the  ELF/VLF  signal .  In  analogy  with  ( 41 ) >  we  can 
write  the  power  downconversion  efficiency  out  to  a  distance 


P3  ( z ) 

pY<0) 


where 


T  dxfia  (O)l^  -  la  (z)l^  +  iL-  o(z)] 

^  -L  ^  ^ _ I _ ^ 

*^1  ?  ~  2 

J  dxla^ (0) I 


~  2 

P  =  ai  la  .1 
j  J  j 


(43a) 


(43fa) 


r(z)  =  ^  IKI  J  dz  Re  (a^j  a^  a^ ) 


(43c) 


Therefore#  in  the  collisional  limit#  the  amount  of  ELF  Power 
one  can  generate  may  exceed  the  collisionless  maximum  Uy/Ui 
set  by  the  Manley-Rowe  relations#  even  when  total  pump 
depletion  does  not  take  place. 


Of  great  practical  interest  to  the  FLF  generation  is 
the  pump  requirement  for  absolute  instability.  From  the 
definition  of  Ic*  we  can  solve  for  the  incident  oumo  field 


'  "  TJT 


;^^2  ^^3  _ ^ 

i^^2x  ^''3x  lAI**  £ 


where  we  have  used  (25)  and  (27)  for  u^#  03#  k>  and  a  in 

(42c).  Recognizing 


=  -2k 


28 


IV  STEADY-STATE  INTERACTION-ELF  GENERATION 


Relevant  to  the  ELF  seneration  scheme  is  the  steady 
state  interaction  siven  by  (26).  With  the  time  variable 
replaced  by  the  initial  value  problem  we  hitherto 
considered  turns  into  a  boundary  value  problem.  The 
absolute  instability  found  previously*  which  corresponds  to 
an  oscillator  in  time*  now  describes  an  amplifier  in  space 
for  the  low-f reauency  wave.  In  the  nonlinear  saturation 
stase*  the  pump  depletes  as  it  propasates.  Due  to  the  fact 
that  (24)  and  (26)  are  equivalent  systems*  it  is  straisht- 
forward  to  see  that  the  conditions  for  absolute  instability 
in  the  steady  state  case  are 


u^U3  <  0 

and 


(42a) 


(42b) 


(42c) 


Since  Ujz’  ^3z  ^  coincides  with  (28a)  while  (41b) 
can  be  shown  to  be  identical  to  (37)*  i.e.  ic  =  Ic*  "ri® 
replaces  the  initial  condition  a^  (t=0)  in  by  the 
boundary  condition  (z-0).  It  follows  that  the 
conclusions  of  the  nonlinear  analysis  on  the  time  dependent 


l(0)  -0,^ 


I ('*1 


(0) l^-la, (t) l^+  —  al 

1  ^3  J 


J  dx  la^(0; 


It  is  easy  to  see  that  in  the  collisionless  limit*  the 
maximum  downconversion  efficiency  one  can  achieve  is  u^/Ui 
if  the  oump  is  completely  depleted*  i.e.»  ajItJ  =  0.  In  the 
example  we  consider*  a^Ct)  ^  O  so  that  the  downconversion 
efficiency  is  less  than  The  collisionless  maximum 
Wa/Wi*  a  consequence  of  the  Manley-Rowe  relations*  can  be 
exceeded  in  the  collisional  interaction  when  the  collisional 
heatins  term  dominates  the  pump  depletion  term*  lai(t)l*. 
Doe  to  the  larse  factor  v/u^,  this  may  occur  when  the  three 
waves  overlap  even  over  a  short  duration.  As  Fisure  4 
shows*  the  downconversion  efficiency  indeed  exceeds  u^/Ui 
without  total  depletion  of  the  pump.  Therefore*  we  have  the 
now  interestins  result  that  in  a  collisional  interaction, 
the  Manley-Rowe  relations  no  loneer  set  the  limit  for  the 
downconversion  efficiency.  The  determinine  factor  in  this 
case  is  the  cumulative  collisional  heatine  over  time.  The 
effect  is  most  sianificant  when  the  convection  of  the  decay 
waves  out  of  the  pump  is  slow  so  that  a  has  the  laraest 
intearated  value  possible.  Under  such  conditions*  tlie 
downconversion  efficiency  exceeds  the  1 1  isionless  maximum 
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the  Manley-Rowe  relations  predict 


the  col 1 i sionless 


interaction  yields  a  value  equal  to  the  frequency  ratio 

=  0.01. 

The  present  numerical  results  can  be  shown  to  be 
consistent  with  the  action  relations  we  previously  derived. 
Subtractine  <38a)  from  <38b)  and  intearatine  in  time>  we 
obtain  (Tj  =  0) 


dx  ( I a^ ( t ) I 


2 


Ia2<t) 1^) 


V 


dx  a 


<40a) 


where 

8  * 

a  =  ^  IKI  J  dt  Re  (a^^  a^  a^ )  (41b) 

o 

For  the  collisionless  interaction  (v  =  0).  (40a)  simply 

states  that  the  transfer  of  action  to  aj  and  33  are 

identical:  hence  the  appearance  of  symmetric  pulses  in 

Fieure  2.  It  then  follows  directly  that  the  enerey  ratio  in 

Fieure  4  should  be  equal  to  Uj/U2>  which  it  is.  In  the 

collisional  case>  (40a)  implies  a  net  sain  in  action- 
transfer  to  33 >  and  we  observe  numerical  evidence  to  this  in 
the  emersence  of  nonsymmetric  pulses  as  well  as  in  the 

increase  in  3(t)/  ^(t).  An  expression  for  the  down- 

conversion  efficiency  may  be  obtained  from  (38b)  by  an 
intesration  in  time  (Fj  =  0) 


different  action  transfer  picture  than  that  described  by  the 
Manley-Rowe  relations.  Me  have  plotted  in  Fieure  3,  the 
percentaee  chanee  of  ^13  time.  Mhile  there  is 
a  moderate  decrease  in  Cj^j'  increase  in  0^3  is  dramatic. 
The  preferential  transfer  of  action  to  33  is  understandable 
in  view  of  the  larse  factor  v/u^  in  the  low  freauency 
equation  as  compared  to  the  small  correction  v/u^  in  the 
equation  for  32-  Physically*  the  enhancement  is  due  to 
additional  (J3  photons  that  are  beine  eenerated  at  a  larse 
rate  via  the  temperature  perturbation  as  a  result  of 
electron  heatins  by  the  h ish-f reouency  waves.  This  is 
achieved  at  the  expense  of  the  total  enersy  7  which  suffers 
a  net  loss  to  the  medium  as  evident  by  the  lona-time 
decrease  shown  in  Fieure  3.  In  comparison*  the 
correspondina  collisionless  interaction  displays  no 
noticeable  chanse  in  any  of  the  three  quantities  considered. 

Of  particular  interest  to  the  ELF/VLF  aeneration  is  the 
efficiency  with  which  one  can  downconvert  the  pump  power 
into  the  low-frequency  wave  &>3 .  Toward  this  end.  we  have 
calculated  the  downconversion  efficiency  3(t)/  ^(0)  for 
both  the  collisional  and  collisionless  interaction.  Fieure 
4  shows  an  order-of -maamtude  (~20)  enhancement  of  the 
former  over  the  latter.  Also  shown  is  the  eneray  ratio 
3(t)/  2^^)*  Aaain  we  find  that  the  collisonal  interaction 
assumes  a  much  laraer  value  ("O.S)  at  pump  depletion.  As 


n.  =0  J  =  1,2,2  (A9) 

whereas  the  nonlinear  temperature  perturbation  is  due  only 

to  Ohmic  heatins* 
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Since  the  thermal  pressure  does  not  affect  the  electron 


velocities  at  hiah  frequencies  =  0)  w®  find 

(Al)  the  nonlinearly  perturbed  velocities 


us  ina 
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where 
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(A12) 


To  arrive  at  (A.L1)>  we  have  combined  the  terms  on  the  riaht 

hand  side  of  (A.l)  with  the  electric  field  term  on  the  left 

hand  side  and  defined  the  nonlinear  fields  which 

contain^  in  addition  to  the  Ohmic  heatina 
ponderomot i VO  contributions  within  the  square  brackets. 

In  a  similar  fashon.  the  nonlinearly  perturbed 

velocity  at  low  frequency  is  aiven  by 
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(A13) 


where 
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Here>  the  ponderomotive  field  is  entirely 
drift  of  electrons.  For  our  oaraineters 
linear  thermal  pressure  n^Tj  turns  out  to 
has  been  ienored. 


+  (A14) 

due  to  the  E  X  B 
of  interest,  the 
be  neelisible  and 


The  nonlinear  current  densities  that  drive  the  hish 
frequency  waves  u)i,  U2  are  due  to  the  nonlinearly  perturbed 
velocities  ^1,2  ^^d  to  the  beatins  of  the  linear  velocities 
loM-requency  density  perturbation 


^  =  -e  (n„v~'-  * 


{Aje>) 
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where  the  asterisk  superscript  denotes  the  complex 
conjusate.  We  keep  only  the  second  of  the  two  contributions 
to  J^l.2  since  it  can  be  shown  that  its  effect  is  sreater 
by  the  ratio  of  h ieh-f reouency  to  low-frequency  phase 
velocities.  Expressins  ^^1,2  explicitly  in  terms  of 

the  electric  fields  via  equations  (A. 4).  (A. 6)  and  (A. 7). 

(A. 15)  become 
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The  loM-f requency  nonlinear  current  density*  however 


is  caused  solely  by  the  nonlinearly  perturbed  velocity*  i.i 


=  -®"o''3 


( A17) 


because  all  hish  frequency  density  perturbations  vanish. 
The  two  nonlinear  fields  that  drive  ♦*'*^3  are  evaluated 
separately.  Substitutins  (A. 10)  for  we  have  for  the 

thermal  field 
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3  3  ®  <^3  i 
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where  Av  has  been  neelected  when  compared  to  (^3 .  If  one 
replaces  by  (A. 4)*  applies  the  frequency  matching 

condition  and  assumes  ~  Ciijp  eJ  may  be  further  reduced  to 
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Similar  treatment  on  the  ponderomot i ve  field 
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where  we  have  related  8^,2  ^i*2  via  Faraday's  law  and 
have  made  simplifications  using  vector  identities.  Dropping 
the  contribution  proportional  to  because  it  is  smaller 
compared  to  ki*  we  finally  obtain 


Fis.  la.— Real  part  of  i?ic  (tjLc)  versus  L/ic  (L/Lc)  for 
the  collisional  (collisionless)  interaction.  The  dotted 
lines  refer  to  the  collisionless  case. 

Fie.  lb.  Imanainary  part  of  tjIq  (1)Lc)  versus  L/Iq  (L/Lc> 
for  the  collisional  (collisionless)  interaction.  The  dotted 
line  at  zero  is  for  the  collisionless  case. 

Fia.  2.  Maves  profiles  at  b=25.0  for  the  collisional 
(solid  lines)  abnd  collisionless  (dotted  lines)  inter¬ 
actions  . 

Fia.  3.  Percentaae  chanae  in  the  sums  of  the  actions 
^12'  ^13  total  eneray  j  versus  time.  The  dotted 
line  indicates  zero  chanae  for  the  collisionless  inter¬ 
action. 

Fie.  A.  The  enerae  ratio  e3(t)/€2(t)  and  the  down- 
conversion  efficiency  €3(t)/€i(0)  versus  time  for  the  colli- 
sional  (solid  lines)  and  coll isionless  (dotted  lines)  inter¬ 
actions. 

Fia.  5.  Geometry  of  three-wave  interaction  in  the  lower 


ionosphere. 


Fis.  6.  “Schematic  of  practical  ELF  seneration  scheme  by 
erazins  incidence  of  two  HF  heaters  in  the  ionosphere. 
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dominance  of  the  EW  magnetic  conqponent.  First,  the  absence 
of  a  vertical  field  component  implies  that  the  nonlinear 
current  flows  horizontally  in  the  ionosphere  in  a  rather 
thin  sheet.  From  Ref.  4  we  note  that  ~  e’  •  ISa  where 

€3  is  the  low  frequency  dielectric  tensor.  Since  the 
dominant  element  of  is  the  off  diagonal  (e^xy) 

current  flows  on  a  plane  perpendicular  to  ka  and  for 
wavepacket  interactions  in  the  horizontal  direction  will 
produce  no  vertical  field.  The  predominance  of  the  EW 
component  over  the  NS  is  consistent  with  the  interaction 
occurring  predominantly  in  NS  direction  which  is  the 
direction  of  the  magnetic  meridian  plane.  Notice  that  in 
the  F-region,  where  the  field  aligned  Pederson  conductivity 
is  dominant,  most  of  the  currents  are  forced  to  flow  in  the 
N-S  direction.  If  the  interactions  would  have  occurred  in 
the  E-region,  where  the  Hall  conductivity  prevails,  we  would 
have  expected  a  different  polarization  for  the  ULF  signal. 
Finally,  the  unsteady  nature  of  the  signals  could  be  either 
due  to:  (i)  We  were  operating  very  near  threshold,  which  was 
attained  only  for  certain  ionospheric  conditions  or  (ii) 
Nonlinear  self  modulation  or  dissipation  effects  forced  the 
system  to  behave  like  a  quenched  oscillator.  Which  is  the 
case  will  require  further  experimentation  amd  better 
diagnostics . 

The  concept  tested  on  the  second  set  of  experiments  was 
similar  to  the  spontaneous  magnetic  field  generation  in 
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where  R  is  the  beam  radius  in  the  F  region  and  17  is  given  by 


w. 


w. 


(3a) 


in  the  collisionless  region  (v  <<  W3)  and 


(3b) 


in  the  semicollisional  regime  (v  >  W3).  For  our  parameters 
Vc  =  10"®  and  Vs  -  3x10”®.  Eq.  (2)  with  Vc  =  lO”*  gives  .32 
my  which  is  smaller  than  the  observed.  For  Vs  =  3x10”®,  we 
find  B3  1.7  my.  On  the  basis  of  the  above  admittedly 
crude  estimates,  we  can  conclude  that  the  signals  produced 
exceed  the  upper  maximum  expected  from  Manley  Rowe  relations 
and  are  in  line  with  the  predictions  of  the  dissipative 
three  way  coupling^'®. 

The  absence  of  informations  concerning  the  precise 
coupling  geometry  and  the  direction  of  the  density 
gradients,  makes  difficult  any  quantitative  comparison  with 
observations  relying  on  the  flow  pattern  of  the  nonlinear 
currents.  We  can  only  comment  on  the  implications  of  the 
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virtiere  (!>«»  are  the  plasma  and  cyclotron  frequency,  L  Is 


the  pump  width  and  A  =  1  +  i  4/3  v/u^  where  v  is  the 
dominant  electron  collision  frequency.  For  W3  <<  v,  and 
k2X/IC3X  *•  .5  we  find 
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For  the  case  under  consideration  the  interaction  occurs  in 
the  P-region  rather  than  the  E-region  considered  in  Refs. 
2,3.  Therefore,  v  corresponds  to  the  electron-ion  collision 
frequency,  which  we  take  as  10^,  Oe/We  •  1/3  and  L  «•  60km 
of  the  order  of  the  density  gradient.  We  also  take  the 
ratio  Wi/Ue  •*  1/sind  *  2  by  using  an  upper  limit  of  ©  ■•  30", 
the  HP  beam  cone.  With  these  values  we  find  Ethr  2-3  V/m. 
For  the  Arecibo  transmitter  with  400  kw  power  and  22-25  db 
antenna  gain  the  field  in  the  F-peak  region  will  be  of  the 
order  of  Ep  =  .3  a  V/m  where  a  is  the  upswelling  factor  in 
the  vicinity  of  the  reflection  area  which  is  of  the  order 
5-10.  Therefore  the  Arecibo  transmitter  fields  in  the 
P-peak  are  of  the  order  required  for  stimulated  beat  wave 
excitation  for  €3  <  5-10  Hz.  On  the  basis  of  achieving 
maximum  theoretical  efficiency*  the  low  frequency  magnetic 
field  will  be  given  by 


four  transmitters  were  turned  on/off  at  equivalent  ULF  rates 
of  5.0  and  6.25  Hz.  There  was  good  coincidence  of  the 
transmitter  operation  and  ULF  detection.  An  example  of  such 
generation  for  5.0Hz  signal  is  shown  in  Fig.  3.  The 
duration  of  this  experiment  was  too  short  to  compile  any 
meaningful  statistics. 

The  limited  number  of  experimental  results  as  well  as 
the  absence  of  diagnostics  of  the  modifications  occurring  in 
the  ionosphere  preclude  at  present  anything  but  the  most 
crude  comparison  of  the  results  with  the  theoretical 
predictions.  For  the  beat  excitation  the  theoretical 
analysis  considered  the  geometry  shown  in  Fig.  4.  Two  HF 
pumps  (oil,  iSi)>  ((^2*  iSa)  were  interacting  to  generate  a  low 
frequency  wave  (0)3,  153)  satisfying  the  conditions  (J1-W2  = 
(J3 .  The  analysis  performed  on  the  basis  of  wave-packet 
interactions  predicts  stimulated  beat  wave  ecitation  leading 
to  pump  depletion  for  HF  threshold  electric  field  in  the 
interaction  region  given  by^ 
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and  the  Schumann  resonance  near  8  Hz.  Fig.  l  shows  an 
example  of  the  data  for  HP  operation  at  5.1  MHz  and  (5.1  MHz 
4-  5  or  3  Hz).  The  HF  transmitters  were  operating  at  a 
difference  of  5  Hz  during  16.30-17.30  AST.  The  frequency 
difference  was  changed  to  3  Hz  at  17.30  AST  and  then  back  to 
5  HZ  at  18.00  AST.  These  data  are  shown  as  power  spectra 
over  the  frequency  0-10  Hz  In  Fig.  1.  The  peaks  correspond 
to  absolute  values  of  150-350  tiy/^z  and  a  signal  to  noise 
ratio  of  about  one.  It  is  worth  noting  that  there  were  no 
communications  on  Mona  and  the  operators  identified  the  3 
and  5  Hz  signals  without  any  prior  knowledge  except  for  the 
prearranged  schedule  that  changes  will  be  made  every  hour  or 
half  hour  and  the  beat  frequency  will  be  between  0-10  Hz. 
All  three  components  of  the  field  were  monitored  during 
these  experiments.  Significantly  though,  the  signals  were 
received  only  in  the  coil  aligned  in  the  East/West 
direction,  while  there  was  no  detectable  signal  in  either 
the  North/South  or  vertical  component  of  the  magnetic  field. 
The  coherent  nature  of  the  signal  can  be  used  to  improve  the 
S/N  ratio  by  using  narrower  filter.  A  significant  finding 
was  that  the  signal  had  considerable  fading.  The  fading 
determined  the  narrowest  filtering.  Using  a  5mH  filter  the 
S/N  ratio  increased  in  occasional  bursts  as  much  as  50  with 
signal  strength  of  the  order  of  1  my //Hz.  This  is  shown  in 
Fig.  2  for  the  3  and  5  Hz  component. 

In  the  second  type  of  experiment  (Feb.  12,  1985),  all 


The  experiments  were  performed  during  the  period  Jan. 
31-Feb.  15*  1985 »  using  the  Arecibo  Observatory  facilities. 
The  HF  transmitters  were  operated  at  5.1  MHz  and  3.1  MHz  at 
a  power  of  4x100  kw  which  was  fed  into  an  antenna  of  22-25 
dB  gain.  The  receiving  system  was  built  around  sensors 
borrowed  from  the  Engineering  Geoscience  Group,  Dept,  of 
Material  Sciences,  University  of  California,  Berkeley.  The 
low  noise  induction  coils  used^  had  a  flat  response  for 
magnetic  field  in  the  frequency  band  .001-500  Hz.  The  coils 
were  portable  and  had  a  sensitivity  with  the  built  in 
amplifier  of  about  500  mV/v.  Early  attempts  to  detect  the 
ULF  waves  at  Arecibo  and  at  Los  Canos  (7km  away  from  the 
heater)  did  not  produce  any  signals.  These  sites  had  large 
background  noise  levels  due  to  household  and  industrial 
circuits.  This  prompted  us  to  move  the  portable  detectors 
to  the  uninhabited  island  of  Mona,  about  150  km  west  of  the 
heater.  The  background  noise  at  Mona  was  6-10  db  below  the 
Puerto  Rico  stations.  We  report  below  the  results  of  the 
experiments  recorded  at  Mona. 

Two  types  of  experiments  were  performed.  In  the  first 
one,  the  four  HF  transmitters  were  split  in  two  pairs  and 
the  desired  ULF  frequency  difference  was  introduced  between 
each  pair.  The  ULF  frequencies  generated  were  3,5  and  6.25 
Hz.  This  ULF  frequency  range  was  selected  by  the  fact  that 
they  are  nestled  in  a  relatively  low  background  region 
between  the  sub  ULF  geomagnetic  fluctuations  below  1-2  Hz 
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Nonlinear  excitation  of  low  frequency  waves  by  the 
interaction  of  two  high  frequency  (HF)  waves  or  an  amplitude 
modulated  HF  wave  has,  in  addition  to  its  intrinsic 
scientific  interest,  many  practical  applications.  it  has 
been  suggested  recently^”®  that  extremely  efficient 
downconversion  from  HF  to  ULF  can  be  achieved  in  the  so 
called  semicollisional  regime,  defined  as  the  region  where 
the  frequency  0)3  of  the  ULF  wave  is  smaller  than  the 
electron-ion  or  electron  neutral  frequency  v  (i.e.  v  >>  W3). 
Two  techniques  were  proposed  in  achieving  this.  The  first 
relies  on  nonlinear  coupling  of  two  HF  waves  ((Ji,lc3), 
((•>2*^2)*  satisfying  the  phase  matching  conditions 
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where  ^  is  the  wavenumber  of  the  ULF  wave.  The  process  has 
a  power  threshold  decreasing  with  W3  and  for  U3  <<  v 
results  in  efficiency  exceeding  that  expected  on  the  basis 
of  the  Manley-Rowe  relations.  The  second  relies  on  the 
spontaneous  magnetic  field  generation*'®  by  coupling  the 
modulated  hot  spot  temperature  gradient  (vT)  created  by  the 
HF  heater  to  the  ambient  ionospheric  density  gradient  (vn) . 
We  report  in  this  letter  the  first  experimental  verification 
of  these  nonlinear  processes. 
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ABSTHACT 


We  report  the  first  experimental  evidence  for 
stimulated  beat  excitation  of  ULF  waves  by  two  HP  pun^s  in 
the  ionosphere,  with  conversion  efficiency  exceeding  the  one 
predicted  by  Manley  Rowe  relations.  Both  the  low  threshold 
for  stimulated  excitation  and  the  high  efficiency  are 
consistent  with  theoretical  predictions  for  helicon 
generation  in  the  semicollisional  regime,  i.e.  when  the  ion 
collision  frequency  is  higher  than  the  ULF  frequency. 
Generation  of  ULF  waves  due  to  spontaneous  generation  of 
magnetic  fields  caused  by  the  coupling  of  the  heater  induced 
temperature  gradient  to  ionospheric  gradient  (vnxvT)  is  also 
reported. 
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laser  produced  plasmas*.  Namely  the  coupling  of  a  density 
gradient  vn  with  a  hot  plasma  region  produces  a  battery  term 
proportional  to  vnxvT,  where  T  is  the  temperature  of  the  hot 
region.  This  acts  as  a  dynamo  producing  a  magnetic  field 
according  to 


(vnxvT) 


As  discussed  in  Ref.  5  the  magnitude  of  the  low  frequency 
magnetic  field  produced  by  an  oscillator  temperature 
gradient  is 
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where  T  is  the  oscillatory  electron  temperature  due  to  the 
on-off  HP  heater  operation,  Ln  is  the  ionospheric  density 
gradient  and  Lr  is  the  horizontal  size  of  the  hot  spot.  For 
our  case  Lm  ■■  50-60km,  Lj  •*  50km  fa  =  5Hz  so  that 
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Namely  mv  fields  can  be  achieved  for  oscillatory  electron 
heating  of  the  order  of  T  -  1000*1^.  Such  heating  is  not 
inconsistent  with  the  one  expected  from  the  available  HP 
power-  Notice  that  Gordon  and  Carlson'’  observed  F  region 
electron  temperature  to  increase  by  about  200*lc  when  the  HP 
power  was  about  4  to  5  times  less  than  the  present 
situation.  As  noted  in  Sef.  5  we  believe  that  this  process 
is  the  cause  of  the  micropulsations  at  .01  and  .02  Hz  with 
amplitude  of  S-lO-y  observed  during  on-off  ‘  heating 
experiments  by  Stubbe  and  Kopka®.  Before  presenting  our 
concluding  remarks  we  should  emphasize  that  both  types  of 
ULF  generation  presented  here  are  different  from  the  low 
frequency  wave  excitation  by  modulation  of  the  polar  or 
equatorial  electrojet*"^® .  The  latter  is  basically  a  linear 
process  which  requires  the  presence  of  ambient  ionospheric 
currents.  The  modulated  electron  heating  modulates  the 
local  ionospheric  conductivity  tensor.  In  the  presence  of 
d.c.  currents,  an  a.c.  current  is  generated  by  conductivity 
modulation  which  radiates  at  the  modulation  frequency. 

In  this  letter  we  have  presented:  (i)  The  first 
experimental  evidence  for  generation  of  ULF  fields  by 
stimulated  beat  excitation  processes.  The  threshold  power 
is  consistent  with  the  thresholds  expected  in  the 
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semicolllsional  regime.  An  iR^ortant  finding  for  nonlinear 
physics  is  the  verification  of  the  prediction^ » 3  that  in 
this  regime  frequency  downconversion  efficiencies  exceeding 
those  predicted  by  Manley-Rowe  can  be  achieved.  (ii) 
Experimental  verification  for  generation  of  ULF  waves  by 
spontaneous  magnetic  field  generation  due  to  the  coupling  of 
hot  spot  oscillatory  temperature  gradients  with  ambient 
ionospheric  density  gradients  (vnxvT) .  Further  studies 
including  scaling  as  well  as  more  precise  measurements  are 
planned  for  the  near  future  and  will  be  reported  elsewhere. 
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FIGURE  CAPTIONS 


Pig.  1  Spectra  of  the  received  signal  In  the  0-10  Hz  band. 

(February  14,  1985).  Receiver  was  located  at  Mona 
Island.  Data  cover  the  period  1630  -  1830  AST.  The 
HF  transmitters  were  operating  at  5.1  MHz  and  with  a 
difference  frequency  Af  of  5Hz  during  1630-1730  AST, 
which  was  changed  to  3Hz  during  1730-1800  AST  and 
changed  back  to  5.0  Hz  during  1800-1830  AST.  The 
magnitude  of  the  5.0Hz  signals  Is  about  160 
and  that  of  the  3.0Hz  signal  Is  about  340  mVHz~^. 

Fig.  2  Spectra  for  February  14,  1985,  during  1737-1740  AST 
and  1821-1824  AST.  The  spectral  bandwidth  is  5  mHz. 
Note  the  large  S/N  ratios  for  both  the  3  and  5  Hz 
signals. 

Pig.  3  ULF  signal  generated  when  the  HF  was  turned  on/off 
at  100  m.s.  rate.  The  receiver  was  located  at  Los 
Canos . 

Fig.  4  Geometry  of  three-wave  interaction  in  the 


ionosphere . 
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Appendix  IV 


"Helicons  in  the  Lower  Ionosphere" 


C.R.  Menyuk  and  K.  Papadopoulos 


HELICONS  IN  THE  LOWER  IONOSPHERE 
by  C.R.  Menyuk  and  K.  Papadopoulos 


Helicons,  waves  which  obey  the  dispersion  relation 
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vfhere  is  the  electron  cyclotron  frequency  amd  We  is  the 
electron  plassui  frequency,  have  long  been  known  to  appear  in 
solid-state  plasmas  for  all  the  way  down  to  w'Oi 
This  wave  appears  rather  than  the  more  familiar  AlfvAn  wave 
because  in  solids  the  ions  are  frozen  into  place  in  the 
lattice,  as  a  consequence  of  which,  the  plasma  current  is 
carried  primarily  by  electrons,  instead  of  by  ions  as  is  the 
case  with  the  usual  AlfvSn  wave. 

In  the  O  and  lower  E  region  of  the  ionosphere  where  z  ^ 
110  km,  one  finds  that  0^,  i.e.  the  ion-neutral 
collision  frequency  is  much  greater  than  the  ion  cyclotron 
frequency. 2  rt  follows  that  the  ions  are  viscously  frozen 
and  only  electrons  respond  to  the  wave  even  when  wzni,  so 
that  helicons  propagate  in  this  region  of  the  ionosphere, 
not  AlfvAn  waves. 
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Ntt  begin  by  recalling  the  electron  and  Ion  equations 
of  motion, 

T.  *  r<«7T»;T  If  *  r.  *  V=1  *  ®' 

(2) 

Yi  -  [f  *  Yi  *  V'l  *  “• 

Here,  m  and  M  are  the  electron  and  ion  masses  and  e  is  the 
ion  charge.  (The  ions,  primarily  NO^  and  0$  in  the  lower 
ionosphere,  are  all  singly-charged).  The  quantities  and 
Vi  are  the  electron-neutral  and  ion-neutral  collision 
frequencies.  Prom  this  equation,  we  find 
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nrriting  i  no«(^i-%)r  %^are  no  la  tha  ion  danaity,  and 
sttbatttuting  into  tha  flald  aquation 


J  ►  "  E 

*  ■ 


(A) 


Mhara  wm  ara  aaauming  for  aimpllcity  that  tha  wave  is 
propagating  along  tha  magnatic  fiald  linas,  wa  now  find 
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Equation  (5)  laada  to  tha  diparsion  relation 
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Vm  now  roduce  Eq. (6)  using  approximations  appropriate 
for  the  region  of  the  ionosphere  vfhich  we  are  considering. 
To  begin  with, 

2  2  2  2  2  2 
wj  '  li.>+iv^l'‘  w“  ,  (7) 

from  %ifhich  wm  conclude 


wJ(<,H>iv^)/«  I 

n^(n^±tt±ie^)  I 


(8) 


%diere  %ra  have  used  the  relation  /  n  -  •  /  O. .  Dropping 
the  displacement  current  contribution  since  we  are 
interested  in  low  frequency  waves  and  using  the  relation 


Vi  Oi  '  w  , 
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we  conclude  through  first  order  in  w/vi  and  n/v^, 
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vdilch  leads  to  Bq.(l)  if  we  neglect  the  small  dasplng 
contr Ibution. 
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This  mod*  appears  ideally  suited  for  BLF  (extremely  low 
frequency)  communication  applications*'*  since  the  lower  in 
the  ionosphere  one  goes,  the  lower  the  t#ave  attenuation 

becomes. 

W*  thank  Or.  K.  Ko  and  Dr.  A.T.  Drobot  for  useful 
discussions.  This  %K>rk  has  been  supported  by  O.N.R. 
contract  •  N00014-84-C-03S37 
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